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SUMMABY 


The  effects  of  conduction  and  viscosity  on  the  stability  of  laminar 
flame  are  examined.  If  <£  denotes  the  ratio  of  the  wave  length  of  a  dis¬ 
turbance  to  the  flame  width  and  o<  is  the  ratio  of  the  \iltimate  temperature 
of  the  burned  gas  to  the  initial  temperature  of  the  fresh  mixture,  the  flame 
is  found  to  be  stable  if 


>L_ 

R.Pr 


where  Re  is  the  Reynolds  number  of  the  flame  based  on  the  flame  width  and 


Pr  denotes  the  Prandtl  number  of  the  mixture.  It  is  further  shown  that 
the  stabilization  is  achieved  primarily  throu^  the  effect  of  heat  con¬ 
duction  on  the  flame  speed  rather  than  the  influence  of  viscosity. 


Effects  of  Conduction  and  Viscosity  on  the  Stability  of  Flame 


1.  Introduction  We  consider  a  mixture  of  two  gases  reacting  with  one  another 
in  accordance  with  the  formula 

Let  Y*  t  respectively  the  molecular  wel£^t  of  the  species  o<  , 

the  mass-fraction  of  the  species  in  the  mixture  and  the  mass-fraction  of  the 
species  at  upstream  infinity.  The  composition  of  the  mixture  can  be  con¬ 
veniently  expressed  by  a  variable  \  defined  by 

Y  = (1.2) 


We  consider  a  low  speed  one-dimensional  flow  with  a  rapid  change  of  tempera¬ 
ture  occurring  in  a  narrow  zone  (width  L)  in  which  the  chemical  reaction  (l.l) 
takes  place.  The  low  speed  condition  is  expressed  by  saying  that  the  Mach 
number  of  the  flame  satisfies 

^  «  '  (1.3) 


We  shall  study  the  stability  of  such  a  flame  by  considering  a  small  perturbation. 

This  allows  us  to  linearize  the  problem.  Let  the  x-axls  be  directed  along  the 

flow,  y  and  z-axes  parallel  to  the  flame.  We  assume  the  solution  is  independent 

1.  V.V  +  wt 

of  z  and  that  its  dependence  on  y  and  time  t  are  of  the  form:  e  ' 


where  k'  is  known  and  represents  the  wave  length  of  the  perturbation,  while 
w  is  unknown.  If  the  real  part  of  w  is  positive,  the  flame  is  stable;  if  it 
is  negative,  it  is- unstable. 

The  flame  profile  is  defined  in  terms  of  the  temperature  rather  than 
the  position  coordinate  x.  It  will  be  preferable  to  use  T  instead  of  x  as  one 


of  the  independent  variables  inside  the  flame. 

2.  flame  Profile  Suppose  that  the  enthalpy  of  the  o<  species  is  given  by 

(2.1) 

where  Cp  is  the  specific  heat  at  constant  pressure,  assumed  to  be  equal  for 
both  species,  and  H*,,  denotes  the  enthalpy  of  formation.  Let  us  define  A  by 

A=  rhifliCHoi"  H*t]  (2*2) 

If  the  Lewis  number  of  the  flame  is  equal  to  unity,  the  flame  profile,  expressed 
as  a  function  of  T,  can  easily  be  shown  to  be  described  by 

”  poo* 

u.=u..l 

when  P  /  P  ;  denote  respectively  the  density,  pressure  and  velocity  of 
the  misture;  the  subscript  "o"  signifies  that  these  quantities  are  associated 
with  the  unperturbed  profile  while  subscript  "oo  "  signifies  the  value  of  the 
quantity  (to  which  it  is  appended)  at  upstream  infinity.  The  flame  profile 
is  not  completely  defined  until  the  position  coordinate  x  is  specified  as  a 
function  of  T.  This  last  quantity,  written  as  X,(T),  depends  on  the  produc¬ 
tion  terms  in  the  flame.  Instead  of  specifying  these  temxs,  it  is  more  con¬ 
venient  to  assume  a  suitable  form  x  such  as  that  shown  in  Fig.  1.  It  is 
clear  from  this  figure  that  the  flame  width  L  cannot  be  defined  exactly. 

However,  its  order  of  magnitude  is  known. 

3. -  Remarks  on  the  Perturbed  Flow  Consider  a  point  X*(T)  on  the  constant 


temperature  surface  T  in  the  flame.  When  the  flame  is  perturbed,  this  point 


moves  to  the  point  Xo(T)i-  X  .  When  the  perturbation  is  small. 


(3.1) 


If  D(T)  is  the  displacement  of  the  constant  temperature  surface  whose 
temperature  is  T,  it  is  clear  that 


X  '^D  ;  L. 


(3.2) 


When  the  flame  displacement  is  much  less  than  the  flame  thickness. 


(3.2) 


it  is  easy  to  see  that  D  is  independent  of  T.  For,  the  flame  thickness  is 


related  to  the  mean  free  path  T  by 
LcJ  ^ 

n 


(3.3) 


where  M  denotes  the  Mach  number  of  the  flame.  Since  the  fluid  dynamic  equations 
apply  only  within  the  continuum  limits,  all  length  scales  must  be  much  greater 
than  .  In.  particular. 


P»»I 


A  comparison  of  (3.2)  and  (3.^)  shows  that  D  is  independent  of  T  so  that  inside 
the  flamer*' 


It  is  possible  to  show  that  for 


P'.'L  Rrwl  !  Re  I 


(3.5) 


(3.6) 


where  Pr  denotes  Prandtl  number  and  Re  denotes  Reynolds  number  based  on  the 
flame  thickness,  viscosity  and  heat  conduction  effects  are  negligible  and  the 
flow  is  unstable.  For  the  viscous  effect  to  influence  the  flame  stability,  we 


*If  D~L  or  greater,  (3.1)  shows  that  (3.5)  is  still  true. 


11 


must  have 

El  f<j  I 

kL 

and,  for  the  conduction  to  affect  the  flame  stability, 

EiB.  M  \ 

V  L 


(3.7) 


(3.8) 


In  order  that  conduction  and  viscosity  may  affect  the  flame  stability  and  a 


continuum  analysis  be  meaningful,  it  is  sufficient  that 


D^L  ;  Re,  R«i  ;  kL  «  \ 


(3.9) 


It  is  true  the  condition  Re,  Pr  1  is  unrealistic.  However,  we  shall  show 
later  that  the  more  significant  case  encountered  in  practice  and  characterized  by 


I)(L  ;  kL<i  ;  R.,B< 


(3.10) 


can  be  considered  as  a  limiting  case  of  (3*9) »  I't  is,  therefore,  sufficient 

to  focus  our  attention  on  (3*9)  which  is  easier  to  handle  from  the  mathematical 

viewpoint . 

•  > 

•  The  condition,  x  =»D=  a  con8^>ant,  valid  inside  the  flame  naturally  indicates 
the  interest  of  dividing  the  flow  field  into  three  regions: 


1)  Upstream:  ^  ) 


2)  Flame: 


"7^  n;  I 

T.. 


3)  Downstream:  T~^T»f.  ^1 


(3.11) 


where  </,  is  the  ratio  of  the  temperature  at  downstream  Infinity  to  that  at 
upstream  infinity.  Since  the  temperature  profile  is  a  continuous  function 
of  X,  the  definition  (3*H)  for  the  combustion  zone  is  not  well-defined.  On 
the  other  hand,  we  may  define  the  combustion  zone  by: 
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1) 

Upstream:  ^ 

^»k. 

2) 

Flame: 

/w  Jb. 

3) 

Downstream:  ^ 

(3.12) 


The  new  definition  tolerates  a  discontinuity  in  derivatives.  However,  such 
a  profile  is  not  so  well  adapted  for  studying  the  effect  of  viscosity. 

4.  Basic  Equations  The  basic  equations  governing  the  motion  of  a  reactive 
mixture  are  well  known.  (See  for  example,  ref.  1  or  2).  These  equations  may 
be  rewritten  in  terms  of  the  new  variables 

T--t 


inside  the  flame.  When  the  fleune  is  slightly  perturbed,  the  disturbances  are 
governed  by  the  linearized  version  of  these  equations.  The  nondimensional 
form  of  these  linearized  equations  is  given  below.  We  shall  use  the  same 
letter  for  a  physical  quantity,  whether  or  not  it  has  a  dimension.  Table  I 
shows  how  the  various  quantities  are  nondlmenslonallzed.  Let  us 

TABLE  I 


Physical  Quantity 

Notation 

Nondimensionalized  by 

Unperturbed  velocity 

u. 

u*. 

Perturbed  velocity 

Position  coordinates 

X,  X, 

L 

Position  coordinate 

n 

k- 

Time 

X 

Composition  variable 

CrT.«/A 

I. 


-7- 


Physical  Quantity 

Pressure 

Density 

Temperature 

Introduce  the  quantity 


TABLE  I  (Continued) 

Notation  Nondlmenslonallzed  by 


f.f. 


P- 

T.. 


L  =  (kL)’  »i 


(4.2) 


which  Is  the  ratio  of  flame  thickness  to  the  wave  length  of  the  disturbance. 

2 

Neglecting  terms  of  the  order  of  M  In  comparison  with  unity,  we  obtain  the 
following ‘System  of  equations:  (cf.,  ref.  5) 

Continuity  Equation 

Momentum  Equations 


... ' 


where 


^r5«i  ’ 

Wurt“i4T)y  ST  £»iTj 

4.1. '  jx»  _  ^  Jji*  _  _L  ,  J<^*l 

irjT*-  $JDSf^3TJ 


(4.5) 


w  -=  V  f* 

Diffusion  Equation 


(4.6) 


4.^ 

'0-x*- 

where  Lewis 


^  ^  61  -V 

s  number  has  been  taken  to  be  unity.  ' 


(4.7) 
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Energy  Equation 


r 

If  we  subtract  (4,8)  from  (4.7),  we  get  a  second  order  homogeneous  equation 
in  ^  .  The  solution  of  this  equation  satiefylng  the  boundary  conditions 

^  *  0  at  both  the  upstream  and  the  downstream  infinity  is  obviously 

^  =  0  (4.9) 

so  that,  when  Lewis  momber  is  taken  as  one,  the  temperature  and  composition 
variables  play  exactly  the  seune  role. 

can  be 

omitted  from  all  the  equations.  Outside  the  flame,  «n  •  so  that 


Notice  that  Inside  the  flame  terms  containing  and 

>T 

<jr  “ 

as  far  as  stability  is  concerned  the  energy  equation  will  intervene  only  inside 


■the  flame  where  it  assumes  the  form 
'^X 


i  T-’- 


c:  0 


T 


(u.io) 


We  thus- see  that  conduction  will  influence  the  stability  of  flame  if.  f?e  ^ JO 
is  of  the  order  of  unity.  The ’mechanism  of  its  action  is  also  clear,  .In 
Landau's  case  (ref,  3)  the  assumption  of  a  constant  flame  -velocity  is  Justified 
as  ^  1  in  which  case  (4,10)  becomes 

•  (Vu)  - 

Eq.  (4.11)  shows  that  when  S/  is  not  too  large,  conduction  modifies  the  flame 
speed-  in  such  a  way  that  the  latter  depends  on  the  flame  curvature.  Thus  we 
have  a  case  similar  to  that  considered  by  Markstein  (ref,  4) . 

5.  Complete  Sdlution  with  Conduction  In  this  calculation  we  shall  neglect 
the  effect  of  viscosity.  We  shall  limit  ourselves  to  the  case: 


R«  rj  I  ^  » 


I  ^  R-  to  I 


(5.1) 
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Inside  the  flame,  we  have 


(5.2) 


^  aT" 


(5.3) 


0 

iT‘ 


ft(u  -  4.  • _ -  - 


(5.4) 


(5.5) 


This  system  of  equations  can  also  be  obtained  by  an  expansion  procedure  in 
a  power  series  of  t'  .  (See  ref.  5)  Since  f  we  see  that  eq.  (5.5) 

implies  eq.  .(5.2),  tiie  latter,  therefore,  need  not  be  considered  further.  If  . 
we  assume  that  x/%^De  ,  eq.  (5.5)>  applied  to  two  sides  of  the 

flame,  gives  the  two  relations:  • 

u, »  ♦  Vs  ^  ]■  • 

where  we  have  reverted  to  the  use  of  quantities  with  dimensions.  Prom  (5.4) 


LVcs  U,.Dk(«<-i) 

Finally  from  (5.3)# 

Cp 


(5.8) 


(5.9) 


Outside  the  flame  the  equations  governing  the  disturbances  are  given  by: 

^  Vy,  _  p 

(5.10) 


'h*‘  ^  5y 

it  -i  X  -(•  •  5x 

h  t  U..bi.  =  . 

“  f.«>y 
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in  the  region  ahead  of  the  flame  front  and  by 


llii  ^ 

4.  U  •  li  2  -L- 

it  ix  f.«iy 


(5.11) 


in  the  region  occupied  by  the  burned  gas.  Assuming  a  solution  which  depends 
on  y  and  t  through  the  factor  e  and  that  Re  u  ^  0  (the  unstable  . 

case),  the  solution  can  be  written  down  without  difficulties  (see  ref,  3). 

At  X  =  o,  these  solutions  yield  the  values: 

p=  A, 
f»=  Ai 

U,3_ 


kA. 


bx 


(5.12) 


kA, 


f*«U««»k  (•^+0 

t  Vj  =  .6^  ^  i - Ai«< 


where 


^-JO- 


U*«lc 


.(5.13)  . 


We  have  to  determine  the  four  unknowns;  A^,  Ag^  Bg^  D  using  the. four  . 
homogeneous  equations;  (5*7)^  (5.8),  (5*9).  Since  A^,  Ag,  Bg,  D  cannot  be  . 
ze^o  simultaneously,  the  determinant  of  the  coefficients  must  be  zero.  This 
provides  us  with  an  eigenequation  for  XI.  Simplifying  the  determinant,  we 
obtain  the  following  equation; 


(5.1»+) 
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If  RePrJC-^w 


(5.1^)  reduces  to  Landau's  equation 


1-  (5.15) 

which  always  admits  a  solution  whose  real  part  is  positive  (unstable) . 

An  examination  of  (5*1^)  shows  that  it  will  have  two  roots  with  negative  real 
parts  if 

ReR  (5*16) 

^  M-l 

which  gives  us  a  stability  criterion.  The  number  2o</(ot-l)  is  close  to  2  for 
large  oC.  It  is  always  bigger  than  2.  Note  that  (5*16)  agrees  with  the  state¬ 
ment  we  made  regarding  (3*6)* 

6.  Extension  to  the  case  Represented  by  (3«10)  As  mentioned  before,  if  we 
should  limit  the  application  of  (5«l8)  to  the  case  i^^l.  Re,  Pr  «  1> 

the  Importance  of  the  stability  criterion  would  be  mainly  theoretical.  We 
now  attempt  to  extend  the  range  of  validity  of  (5.I8)  by  reducing  eC  .  The 
decrease  of  X  will  have  two  effects: 

1)  New  terms  will  become  important  in  the  equations. 

2)  It  will  alter  the  locations  where  the  boundary  conditions  are  applied. 
Of  course,  an  intrinsic  difficulty  associated  with  a  small  JD  (s.g.,  n  l) 

-1 

is  to  know'  what  happens  to  our  assumption  of  a  fluid  continuum.  Since  D  K  , 
the  displacement  of  the  constant -temperature  surfaces  becomes  smaller  and 
smaller  ns  X— »1.  In  such  a  case,  it  should  have  no  effect  at  all  on  the 
flame  so  that  the  term  inside  the  flame  can  again  be'  neglected.  As  a 

consequence,  the  energy  equation  inside  the  flame  is  still  given  by  (5*5)* 

The  equation  (4.9)  is,  of  course,  still  valid.  We  shall  limit  ourselves  to 
the  study  of  conduction  effect  while  ignoring  the  viscous  effect  so  that  the 
viscous  terms  in  the  momentum  equation  need  not  be  considered  for  the  moment. 

As  Jti  decreases,  terms  of  the  form  JL  U  . . . .  appear  in  the  continuity 

X»  4T  ' 
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and  momentum  equations.  Since  these  terms  invariably  contain  the  factor 
which  is  of  the  order  of  !/•<  ,  these  terms  are  actually  of  the  order  of  l/«Xi 

Now  in  an  actual  flame,  oi.  is  of  the  order  of  10  so  that  unless  cD  becomes 
too  small,  these  terms  can  again  be  neglected.  It  thus  follows  that  the 
equations  valid  inside  the  flame  remain  the  same  as  those  given  in  eqs.  (5*2) 
to  (5.5).  The  equations  governing  the  disturbance  outside  the  fleime  are  of 
course  unchanged  when  viscous  effect  is  suppressed. 

Next  we  examine  the  effect  of  decreasing  the  quantity  X  on  the  locations 
where  Lhe  boundary  conditions  are  applied.  The  solution  valid  outside  the 


flame  contains  terms  of  the  form  e‘ 


and  e 


.  This  solution  is 


joined  to  the  solution  inside  the  flame  at  the  two  borders  of  the  flame. 

Since  X  is  of  the  order  of  L  or  D,  all  the  exponentials  will  tend  to  I  if 
X  1.  As  X  is  decreased,  this  simplification  does  not  apply.  Taking 
the  origin  of  our  coordinate  system  at  the  middle  of  the  flame,  the  boundary 
conditions  must  now  be  applied  at 

,  ^  iky  +  ot 

(6.1) 

tkx 

When  this  is  substituted  into  the  exponentials  e  and  e  •<  ,  the  result 

can  be  written  as  a  product  of  two  terms  corresponding  to  the  two  terms  in 
(6.1).  The  term  containing  D  will,  of  course,  tend  to  1,  since  D  «  k'^. 

The  term  containing  L  does  not  present  any  difficulties  as  it  enters  into 
the  formulas  as  a  coefficient.  The  fact  that  L  is  not  well  defined  does  not 
matter  either,  since  the  coefficients  are  finally  eliminated  and  do  not  enter 
into  the  expression  of  the  stability  criterion.  In  conclusion,  we  see  that  the 


stability  criterion  (5«l8)  is  valid  for  X^l« 


7.  Stable  Mode  It  is  Interesting  to  calculate  the  form  of  the  stable  mode. 


Instead  of  assxjmlng  Re  (rj.)>  0  ,  we  take 

R4n)<  0 


(7.1) 
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The  solution  of  the  system  of  equation  outside  the  flame  can  he  constructed 
as  before.  At  the  boundary,  we  have 


This  equation  shows  that  if 


K*Pr4<»  (7.13) 

there  is  at  least  one  stable  mode.  To  be  more  precise,  there  will  be  only- 

one  stable  mode  if  RePrJt  <  while  two  such  modes  are  found  for 

M-l 

RePr  greater  than  _  but  less  than  unity.  Finally,  if 

^  *  I 


•i-t 


(7.14) 


the  solution  is  still  stable  since  there  is  no  unstable  solution  for  RePrJC/<_lsL 

¥-1 

Here  n.  is  purely  imaginary  and  the  detailed  structure  of  the  disturbance 
at  any  instant  can  be  calculated  by  solving  an  initial  value  problem. 

8.  Effect  of  Viscosity  Here  we  assume  the  conduction  effect  can  be 
neglected.  We  shall  consider  the  case 

(8.1) 

Applying  a  similar  reasoning  as  used  in  the  previous  sections,  we  find  that 
the  equations  valid  inside  the  flame  are 

(8.2) 


A 


\  -  >  r  I  Vv  I  »  1 


^  V 


JT 


(8.3) 


(8.4) 


where  v  is  given  by  (4.6).  Substituting  (4.6)  into  (8.3)  and  meeting  us.^  of 
(8.4),  we  obtain: 

J-  IL -i- -i—  (8.5) 

which  may  be  used  instead  of  (8.3). 

Outside  the  flame  we  must  use  the  linearized  form  of  Navier-Stokes 
equations.  Thus,  instead  of  (5. 10)  and  (5.11),  we  have: 


(8.6) 
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'*rF  '^x  >y  •  Sx*"  >/*■ 


(8.7) 


(8.8) 


where  f.  ,  stand  for  P  ^  ^  U„ 


in  the  region  ahead  of  the  flame,  while 


they  stand,  for  t./«  ’  in  the  region  behind  the  flame.  The  ap¬ 

propriate  solutions,  applicable  to  the  upstream  and  downstream  regions,  of 
the  above  system  can  be  easily  constructed.  To  obtain  a  stability  criterion, 
we  need  only  the  values  of  the  various  flow  variables  at  the  boundaries  of 
the  flame.  At  the  border  adjoining  the  region  occupied  by  the  unburned  gas. 


VI,* » Ai  /  t) 

ivts  A,/j^»U*«(D.+i) 

At  the  border  adjoining  the  region  occupied  by  the  burned  gas, 

i  r  Ai«<  / 
uvjs - 4. 

where  Xi  is  approximately  given  by 

JL- 

•< 

Now  from  (8.2),  (8.U)  and  (8.5),  we  have 
U],  U,  * 


(8.9) 


(8.10) 


(8.11) 


(8.12) 


Substituting  (8.9)#  (8.IO),  (8.II)  into  (8.12),  we  obtain  a  system  of  homo¬ 
geneous  equations  in  A^,  Ag,  Bg  and  D.  Setting  its  determinant  equal  to  zero. 
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we  have 


f\j  0 


(8.13) 


One  of  the  solution  Is 


D-: 


I  * 


(8.1U) 


as  compared  withfl  in  Landau's  case.  Hence,  a  necessary  condition 

for  stability  is 

(8.15) 

A  comparison  of  (5.l6)  and  (8.15)  shows  that  the  flame  is  primarily  stabilized 
by  the  influence  of  heat  conduction  rather  than  the  effect  of  viscosity. 
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